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Abstract 

On the basis of the one-band t-t'-Hubbard model a self-consistent renormalization the- 
ory of magnetic short-range order (SRO) in the paramagnetic phase is presented combining 
the four-field slave-boson functional-integral scheme with the cluster variational method. 
Contrary to previous SRO approaches the SRO is incorporated at the saddle-point and 
pair-approximation levels. A detailed numerical evaluation of the theory is performed at 
zero temperature, where both the hole- and electron-doped cases as well as band-structure 
effects are studied. The ground-state phase diagram shows the suppression of magnetic 
long-range order in favour of a paramagnetic phase with antiferromagnetic SRO in a wide 
doping region. In this phase the uniform static spin susceptibility increases upon doping 
up to the transition to the Pauli paraphase. Comparing the theory with experiments on 
high-Tc cuprates a good agreement is found. 



PACS number(s): 75.10.-h, 71.28.+d, 71.45.-d 



1 Introduction 



In the theory of strongly correlated itinerant electron systems two topics are of continued in- 
terest, namely (i) the localized-itinerant complementarity and (ii) the competition between 
magnetic long-range order (LRO) and short-range order (SRO). In particular, the concept 
of magnetic SRO and its interrelation to itinerant properties was investigated in the context 
of both narrow-band magnetism of transition metals and their compounds [1-6] and of the 
unconventional magnetic behaviour of high-Tc copper oxides [7-12]. In the cuprates, neutron 



scattering [13| and nuclear magnetic resonance experiments [|l^ reveal pronounced antiferro- 
magnetic spin correlations within the Cu02 planes which persist even in the superconducting 
phase. Moreover, measurements of the spin susceptibility xi^, S) in the normal metallic state 
of La2_5Sr5Cu04 (LSCO) [||, |l| show a maximum in the doping dependence as well as (for 
6 ^ 0.21) in the temperature dependence and give further evidence for strong SRO effects at low 
temperatures, where the SRO decreases with increasing doping and temperature. Thus the ex- 
periments on cuprates bring out the importance of electron correlations as compared with the 
situation in traditional band magnetism, and therefore yield a new challenge for a microscopic 
theory of SRO in itinerant electron systems. Such a theory has to provide a self-consistent 
description of strong SRO (in the absence of LRO) down to zero temperature. 

The previous theoretical approaches to the problem of SRO in transition metals and 
cuprates are mostly based on Hubbard-type models [1, 2, 4-6, 8-12, 17]. For example, to 
explain the normal-state susceptibility of LSCO in terms of SRO, the three-band Hubbard 
model was used In connection with the study of SRO some general problems have gained 
a renewed interest, such as the stability of various magnetic LRO phases against paraphases 
with and without SRO as well as the question of phase separation in strong correlation models 
11,0]. 

Prom the methodical point of view, the SRO theories for itinerant systems are preferen- 
tially formulated within functional-integral representations of the one-band Hubbard model 
using the static approximation. Thereby, various Hubbard-Stratonovich two-field methods, 
often combined with the single-site coherent potential approximation (CPA) [j|, |2|, ||, |8| and, 
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in the context of high-Tc's, the scalar four-field slave-boson approach |19, ^, 12] where em- 
ployed. Those theories are designed to describe the formation and ordering of local magnetic 
moments in an itinerant system on a time scale large compared to the electron hopping time. 

In the mode-mode coupling theory of spin fluctuations (SF) by Moriya et al. [|l|, ^ inter- 
polating between the weakly magnetic (local SF in g-space) and local moment (local SF in real 
space) limits, the SRO is reflected in the spatial correlation of thermal SF (with an amplitude 
increasing with temperature) and is appreciable even above the Curie temperature. Starting 
from the local moment limit and using the two- vector-field Hubbard-Stratonovich/CPA ap- 
proach, the SRO is taken into account by an expansion in pairwise terms (within the bilinear 
approximation) around the single-site CPA saddle-point. The lack of self-consistency in this 
approach may be justified when the SRO is weak. Note that, at T = 0, the theory by Moriya 
et al. ||l|, ^ reduces to the Hartree-Fock approximation in both the weak- and strong-coupling 
limits, and the SRO is lost. The neglect of important correlations at zero temperature may be 
the reason for the absence of a maximum in the calculated spin susceptibility. Contrary, the 
phenomenological version of the interpolation theory given in Ref. [||] incorporates the SRO at 
the saddle point, and in the weakly magnetic limit the results of the classical approximation 
to the self-consistent renormalization theory of SF are recovered. 

In the local-band theory of ferromagnetism where the fiuctuations of the local 

magnetizations with a fixed amplitude are considered to be local in g-space, the SRO is de- 
scribed by an expansion around the Stoner saddle point [Q] and turns out to be strong near 
the Curie temperature. As in Ref. Q, the SRO-induced renormalization of the saddle point 
is disregarded. 

The non-self consistency is also inherent in the SRO theories taking the SF to be local in 
real space and employing various cluster methods |^, |9|. For example, in the approaches of 
Refs. Q and resulting in an effective Ising free-energy functional, where the exchange energy 
is evaluated beyond the bilinear approximation used in Ref. and is nearly independent on 
temperature, the SRO is treated by an expansion around the single-site CPA saddle-point and 
by the Bethe-Peierls approximation 



3 



To sum up, all previous microscopic theories based on Hubbard-Stratonovich or slave- 
boson/CPA approaches do not take into account the SRO self-consistently. Just that, however, 
is needed when we are dealing with a strong SRO at low temperatures. 

Besides the functional-integral techniques, there are only a few alternative methods con- 
cerned with the description of SRO [|l^, |ll[. For example, to discuss the doping dependence 
of the zero-temperature susceptibility x(0, (^) in the one-band Hubbard model by a semi- 
phenomenological weak-coupling approach it was argued that the susceptibility of an 
antiferromagnet with only SRO can be approximated by the perpendicular susceptibility of an 
antiferromagnet with LRO. 



In previous work |12] we have briefly sketched an improved (non-CPA) treatment of SRO 
in the paraphase of the two-dimensional (2D) one-band Hubbard model, based on a four-field 
slave-boson (SB) functional-integral approach, by taking into consideration the feed-back 
effect of SRO on the saddle point. In the present paper our self-consistent renormalization 
theory of SRO, being valid also at T = 0, is presented and discussed in more detail. As 
compared with Ref. [|l2| , the theory is improved methodically, and several extensions and new 
results are given. For example, we investigate the influence of the band structure and of both 
electron and hole dopings on the stability of magnetic SRO versus LRO and on the uniform 
static spin susceptibility in the paraphase. 

The basic idea of our scenario for the occurrence of SRO in itinerant electron systems 
may be understood as follows. Measuring the 'localization' of the electron spins by the lo- 
cal (not thermally induced) magnetic moment, with increasing correlation strength, i.e. with 
increasing U/t and decreasing doping, the local moment increases. Concomitantly, local mag- 
netizations develop and an effective exchange interaction will tend to align neighbouring local 
magnetizations antiferromagnetically, even in the absence of LRO. Of course, the occurrence 
of SRO depends on the doping level, and SRO effects become strongest at T = 0. That is, 
if the local-moment behaviour is sufficiently pronounced as compared with the itinerancy, a 
paraphase with antiferromagnetic SRO becomes energetically favourable. Let us emphasize 
that our SRO concept differs from that described, e.g., in Refs. |jl], |2[. 
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The paper is organized as follows. In Sec. 2, the action of the SB functional integral 
for the partition function of the one-band i-f'-Hubbard model on a square lattice (which 
may describe the essential features of low-energy charge and spin excitations in the cuprates 
|l7|) is expressed in terms of fluctuating local magnetizations and internal magnetic fields 
and transformed to an effective generalized Ising model with exchange integrals of arbitrary 
range. The SRO is taken into account by the summation over the Ising spins within the 
cluster variational method (CVM) [^]. In Sec. 3, the resulting effective free-energy functional 
is treated in a non-local saddle-point approximation, where the SRO is incorporated, at an 
equal level of variational approximation, in a fully self-consistent way. Sec. 4 is devoted to 
the calculation of the uniform static spin susceptibility, where the analytical results on the 
influence of SRO are given within the Bethe-Peierls approximation. In Sec. 5, at T = 0, 
numerical results for the local magnetic moment and magnetization amplitude, the ground- 
state phase diagram, and the spin susceptibility as functions of both the hole and electron 
doping level and of the interaction strength are presented. Thereby, the influence of the band 
structure, as parameterized by the nearest- and next -nearest-neighbour transfer integrals, is 
studied. The summary of our work can be found in Sec. 6. 

2 Effective free— energy functional for the Hubbard model 

2.1 Slave— boson functional integral 

In the scalar four-fleld SB representation [|^] the Fock space at each site is enlarged by in- 
troducing the Bose flelds e,, pi^ and di describing projection operators onto empty, singly and 
doubly occupied states, respectively, and the Hubbard model is expressed as 

w = E %4/l/,.^,. + ^ E - /i E ^fU^^ ' (2-1) 

ij(T i i(T 

where tij are the transfer integrals between nearest (— t) and next-nearest (— t') neighbours, h 
denotes the uniform external magnetic held, and 

Zia = (1 - 4di-pLPiayH4Pia+pl-A)(.'^ - - pVaPi-a)~^ ■ (2-2) 
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To exclude unphysical states in the extended Fock space, the local constraints 



cr 

p\aPia + dUi = flfia 



(2.3) 



(2.4) 



have to be fulfilled. Expressing the partition function Z hy a coherent -state functional integral 
over complex Bose {{b^*^}; b = e, pa, d) and pseudofermionic Grassmann ({/,/}) fields, the 
local constraints are enforced by the time-independent Lagrange multipliers A^^^ and Xj^ ■ 
Integrating out the pseudofermions we have 

Z = |[MW][PeW][PpW][dA«][dA(2)]exp{- j^''dr£(r)}, (2.5) 



+d:(9r + A«-EA(? + ^)d,-A« 



Tr In 



G^lir)] (2.6) 



with 



G^jiir) = i-dr + - XrJ + ah) 6,^ - z^^t,^ . (2.7) 
To remove the non-diagonal fluctuations in the transfer term of the inverse Green propaga- 



tor (|2.7D, we employ the Shiba transformation [22 



Im 



(2.i 



where z^j^ = z^^5ij. Then, in the frequency representation of the Lagrangian (2^), under the 
trace operation G^-^(u;) can be replaced exactly by 



G,^-i(^^ + /u) 



(2) 

uj + fi- XlJ + ah 



5ij tij 



(2.9) 



with qi, 



Zin 



By a local {/(l)*^^ gauge transformation the phases of the Bose fields Cj 



and pia are removed (radial gauge), and the Lagrange multipliers become time-dependent 

(2) 

fields. Applying the static approximation for the bosons {p^^ , , A^^ , where the saddle-point 



approximation for x\^^ is used to eliminate the integrals over the Cj fields) and adopting the 
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transformed fields |^| 



cr (7 



(2.10) 
(2.11) 



the partition function is obtained as 



[Vd] [Vd*] [Vn] [Vu] [Dm] [V^] exp {-/3l^({d, d* , n, v, m, S})} , 



(2.12) 



^ = (Ud*di - mvi + miCi) + M dojfiio - f,) Im Trln \-G^'M] ■ (2.13) 



In eqs. ( |2.1[l| ) and ( 2.11 ), and Ui are the bosonic representations of the local magnetization 



and particle number, respectively, defined analogously to its fermionic counterparts m{ = 
J2a^ficrfia ~ J2c7 fiafia' Siucc .^j couplcs to TTij as a magnetic field we denote by 

'internal magnetic field'. 



2.2 Transformation to a generalized Ising model 

To incorporate the SRO beyond the uniform paramagnetic (PM) saddle point, we perform an 
expansion in terms of the local perturbation V^^S^^ = —G^ja + G^-J, where 



„o „ik{Ri-Rj) 



k ka 



is the PM saddle-point propagator with 



(2.14) 



k ' 



(2.15) 



A° = — 0"^°, and the 2D tight-binding band dispersion 



er = — 2t(cos + cos ky) — 4t' cos k^ cos ky 



(2.16) 



By the diagonal part of the Green propagator ( 2.14 ), the perturbation reads 



(2.17) 
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and we can split up the fermionic part of ( |2.13| ) as 
Tr In 



Tr In 
+Tr In 
+Tr In 



- G%)MTjA^) 



(2.18) 



Here, T^^ = V^^ [1 — G^j^Vj^]^ is the scattering matrix, and G"ijy denotes the nondiagonal part 
of ( p.l4| ). The second term on the r.h.s. of ( |2.18| ) yields a single-site fluctuation contribution to 
the functional ( p.l3|) , whereas the third term describes the coupling between the fluctuations 
at all sites and is responsible for SRO effects. 

Treating the fluctuations of the local magnetizations and the internal magnetic fields 
on an equal footing, we express those fields by their amplitude and direction according to 



± . 



(2.19) 



Furthermore, we make the ansatz 



bi —I- with b = m,^,, n, u, d, d* 



(2.20) 



which becomes relevant for the calculation of the uniform static spin susceptibility (Sec. 4); 
for h = 0, we have bg^ = b. In Appendix A the resulting partition function is transformed to 
that of an effective generalized Ising model along the lines indicated by Kakehashi [0. In the 



pair approximation, where all terms of (2.18) involving more than two sites are neglected, we 
obtain 



with 



Z = Y. /[Pda][M;][Pn,][Pi.,][Pm„][Pea]exp{-/3.P'(6,,{sJ)} (2.21) 



^{bo^,{s^}) = •^(6a)+'^'(P"''^)(fea,{s.}) 



(2.22) 
(2.23) 



Here, a = ±1 and the index n denotes the n-th neighbour shell around a given lattice point i, 
where {ij)n sums over the z„ sites j in the n-th shell. The expressions for h and the effective 
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Ising-exchange integrals J„ are given in Appendix A. Let us emphasize that those quantities 
are comphcated functions of all SB fields and have to be determined self-consistently at each 
interaction strength U and hole doping 5 = 1 — n. 

Furthermore, remark that J„ has to be contrasted from the effective Ising-exchange en- 
ergy J occurring in previous Hubbard-Stratonovich/CPA approaches |6|, 24], where J couples 
thermally induced local moments. 



2.3 Cluster variational method 



To perform the Sj-sum in the partition function ( ^.2lD of the effective Ising model ( 2.23 ), we 
restrict ourselves to the square lattice and employ the CVM in the version given by Hijmans and 



de Boer [21|. In this method the calculation of the free energy in the (t(;)-cluster approximation 
starts from the choice of a basic cluster (w) and the series of subclusters (v) generated by 
overlapping two basic clusters in all possible ways, as illustrated in Fig. 1 [(w) = (/), (v) = 
(e), . . . , (a)]. In the same way, the next generation of overlap figures can be obtained from each 



basic 
cluster: 



(f) 





(d) 



(a) 

Fig. 1: Genesis of the series of subclusters (e) to (a) from the basic cluster (/) in the square 
lattice (see text). 

of the subclusters (v). Note that in each generation there occur 'essential' overlap figures [e.g., 
(e) in Fig. 1] which do not appear as overlap figures in the subsequent generations. 



Following Ref. [21|, the free energy of the A^site system in the (i(;)-cluster approximation 
is obtained as 

w 

•^«.(pi"^) = Eyi"^-^^'^nP?), (2.24) 
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where 

t^^Xp^:^) = ny: Ai^vi^) (e^:^ + ^ inp(:A (2.25) 

K€{v) VP/ 

is the free energy of an assembly of independent clusters of type {v). Here, k labels the 
topologically inequivalent spin configurations of the cluster (v) with energy eii^ and the multi- 
plicity xli^ (number of equivalent configurations of type k), where el^^ and aI^^ are calculated 
for an isolated cluster (v). The distribution number pl^^ is the number of (v)-clusters with 
fixed configuration k contained in a given configuration of the A^-site lattice divided by the 
total number of (f )-clusters in the lattice, where the normalization condition 

E Al'Vl'^) = 1 (2.26) 

has to be fulfilled. Of course, the distribution numbers pi"'* for the subcluster (v) are determined 
by the distribution numbers pl"'^ for the choosen basic cluster (w), i.e. 

Plr^=pi'^^b?], V(^)C(^). (2.27) 

Since in the cluster construction of the total free energy in terms of the 'microvariables' 
pi"'* [p|.T'*] described so far the different subclusters (v) are counted several times, in ( p. 24 ) 



(v) 

the geometrical factors yw were introduced to compensate this overcounting and to give the 
correct number of (t;)~clusters contributing to the free energy. In Appendix B the consistency 
relations ( |2.27| ) as well as the linear system of equations for are given in a more explicit 
form for the most simplest basic clusters (a), (6), and (c). 

Finally, in the spirit of the CVM, the free energy has to be minimized with respect to the 
'macrovariables' p^k'^ Since the constraint ( 2.26| ) also holds for the basic cluster (w), we can 



reduce the variational degrees of freedom by one introducing, instead of the Pk"\ a new set of 
independent variables x^. Accordingly, the CVM equations read 

dfJ^%^^^_ (2.28) 

ox, 

It is worth emphasizing that alternatively all the CVM variables contained in J^y^ can be 
expressed by the complete set of M~point (M = 1, . . . , A^*-"'^) correlation functions K belonging 
to the basic cluster (w) [ESl. 
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Applying this method to our spin model ( |2.23| ) we obtain the effective free-energy func- 
tional in the (it;)-CVM approximation (c.f. Appendix B), 



1 Xi 



-N 



where the internal energy part is given in terms of the expectation value 

and of all pair-correlation functions in the basic cluster, 

2 dTyj 



Kj := {siSj)n 



ZnN dJn 



(2.29) 



(2.30) 



(2.31) 



In (2.29) the coupling to higher, e.g., 3-point or 4-point correlation functions results from the 
entropy term. 



3 Saddle point with short— range order 



Now we calculate the free energy per site 



f{n,h,T) 



SP 



where 



S = [3 



!^(6„) + a^(P-'^)(6,, 



(3.1) 



(3.2) 



is the effective bosonic action, adopting the saddle-point and (ty)-cluster variational approxi- 
mations for all Bose fields ha and variational parameters x^, respectively. Taking into account 
the dependences on of h and. (cf. Appendix A), from dS / db(^ — we get 



n,—a'a'(T 



dbn 



. (3.3) 
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This coupled system of self-consistency equations determining the saddle point with SRO can 
be cast to the form 

m, = aJ2anL , 4 = -^g/^^, (3.4) 



'^a — ^ "-ao- 5 — ^ Qaa p. i (3-5) 



Here, n^^ is the pseudofermionic expression for the particle density given by 



/ 



\ y\(l + a'K 



+E 



T 



(3.7) 



where Xal = ^'a ~ ot^^a- ^l^- ( |3.4| ) to ( p.6p , Q^g- is given by ( |3.7| ) with dXal replaced by 
dqaa- Finally, the chemical potential is determined from the number condition 

n = l-5 = Y^p^^KL (3.8) 

OCT 

with 



pL"^ = ^(1 + «^i)- (3.9) 
At the saddle point, we have n^^ = = {n^+aam^) /2, where n^^ is the bosonic expression. 
For vanishing local magnetization, fri^ = 0, we have = and (by ( [All ) and (|A1^ )) J„ = 
so that there is no SRO, and the PM saddle point (with the external field h) is recovered from 



(3.4) to (|3.6| ). Accordingly, in the /i = limit (m^ = m] i.e., the integrals # in (|3.7| ) only 
depend on = acr = it) we obtain two possible paramagnetic phases: (i) the paraphase 
without SRO (PM) and (ii) the paraphase with magnetic SRO (SRO-PM), i.e., 

PM : (s,) = , {s^s^ = ; m = , 

" ^ (3.10) 

SRO-PM : (sj = , {s^Sj) / ; m > , 
where i and j are nearest-neighbour sites. 

Let us stress that, at T = 0, the SRO-PM phase with antiferromagnetic correlations 
((SjSj) < 0; J]^ < 0) must be distinguished from the phase with antiferromagnetic LRO (de- 
noted by AFM) having a finite sublattice magnetization, = ~ P\i ~ ~^b^ which is 
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determined from the A-B saddle-point solution |26] and differs from the SRO amplitude m; 
in particular, we found m 7^ in parameter regions, where rriA = (see below. Fig. 2). 

4 Spin susceptibility 

The uniform static spin susceptibility x = lim/i-»o ^1 where m = J2aPa''fna with m„ = arh^ 



is the averaged magnetization, has to be calculated from 



using the solutions of the saddle-point and cluster variational equations for and pa\xi), 



(^g^^dnia , dp, 



(a)- 



respectively. The first term in ( [4.1D describes the change of the magnetization amplitude with 
the applied magnetic field and gives mainly the 'itinerant' contribution to x- The second term 
describes directional fluctuations of the magnetizations and is called the 'local' contribution 
being finite only in the SRO-PM phase. However, note that the 'itinerant' and 'local' properties 
are interrelated and determine both contributions to the spin susceptibility. 
Now we exploit the symmetry relations in the paraphase '^^j^ 



dh 



= m, = —m and 

h=0 -h 

= ^ = — ^_ . Taking into consideration the full dependence on the external field 

h=0 h h 

h of the SB fields [rhj^, m°, and variational parameters [xj, in the limit /i — > the 
uniform static spin susceptibility can be obtained from the solution of the following linear 
system of equations 

Sx.y^if^ + ^rh^j^ -]y^+h ^^[+ ^]y^+h ^ '^rn°y^h + '^^"j/Ch = ~^hy ) (4-2) 



S^oyorrff^ + S^oyoCh — ~^hy° ) (4-3) 



where Xy = ^ 



h=o ' \+,-] - ( Ox ax ) ^ 



dxdy 
h=0 



with y = 



and y° = m°, 

In the paraphase without SRO, from (|4.3| ) we find 



= ^ + 1 (44) 

l + Ax°o + Bx'l + lB^xf-X"oX°2) ' 



and the suceptibility is 



l + Ax°o + Bx1 + \B^xf-X°oX"2) 
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with 



A = y^ 



and the higher-order Lindhard functions 

1 



Xi 



N 



dE°. 



(4.6) 



(4.7) 



Let us point out that our result for the spin susceptibihty x° agrees with the static and 
uniform hmit of the dynamic spin susceptibihty derived, within the spin-rotation-invariant 



SB scheme [27|, from the Gaussian fluctuation matrix at the PM saddle point pM]. Contrary. 



the uniform static spin susceptibility given in Ref. |29] within the scalar four-field SB approach 
disagrees with the result of Li et al. p8[ | away from half-filling. 

If we include magnetic SRO effects at the Bethe-Peierls level of approximation {i = 1,2; 



cf. Appendix B), from the evaluation of (4^) together with ( [4.4D - (|4.7] ) we finally obtain the 
spin susceptibility as 

, = , (4.8) 

where 



Xi 



Xii 



A 



'[+,■ 



,«4 



(4.9) 



f3fh+r + /3Z\ 



-r5 



+ 1 rSx^fhj^ ~ "^+^^mj^ 



[+,- 



011 



2exp{-2/3Ji} - 1 



(5A 



, (4.10) 
(4.11) 

(4.12) 



and 



A 



n -1 



-s, 



hi. 



(4.13) 
(4.14) 
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^ — ~Shm^ ~ •Sm°m^''^h ~ ^^°m^^h^ (4-15) 

r = hh + hmoml + h^oil. (4.16) 
In this paper we are primiliary interested in the zero-temperature susceptibility. Equa- 



tions to ( |4T6D then lead to 

Xi for J < (antiferromagnetic SRO) 



lim Y 



X° for J = (PM) (4.17) 

— for J > (ferromagnetic SRO). 

It is noteworthy that, in contrast to the theory by Baumgartel et al. we get a finite spin 
susceptibility at T = 0. 

5 Numerical results and discussion 

We have numerically solved the self-consistency equations ( |3.4[) to ( p. 91) at T = 0, where the 
SRO is included within the Bethe-Peierls approximation, and the 2D tight -binding unper- 
turbed density of states corresponding to ( p. 161) is used. To determine the magnetic ground- 
state properties (SRO versus LRO, susceptibility x(0, (5)), all the SB fields ha and CVM vari- 
ables (j;i, X2) are evaluated as functions of U /t and doping. Note that in the tedious numerics 
particular attention has to be paid to the analytical behaviour of the complex logarithm ap- 
pearing in the integrals ( [A.5| ), (A. 15) and in their derivatives. 



The band structure is incorporated via the non-interacting dispersion ( p. 161) , where we 
have chosen (besides t' = 0) the ratio t' /t = —0.16 according to ab-initio parameters for hole- 



doped {5 > 0) LSCO and ARPES data on electron-doped {6 < 0) Nd2-5Ce5Cu04 fsg]. 
For both compounds, realistic values oiU/t lie between 6 < C//t < 10 ||T7|| . 

5.1 Long— range order versus short— range order 

Concerning the LRO ground states of the 2D t-t'-Hubbard model, within our scalar four- 
field SB approach resulting in an effective Ising free-energy functional, only antiferromagnetic 
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(AFM) and ferromagnetic (FM) LRO may be described by the A-B and uniform saddle-point 
solutions, respectively pl| , p6[ . Up to U/t = 10 and for either doping level we do not obtain 
a FM phase. This should be qualitatively contrasted with the results obtained by Hartree- 
Fock approaches p2| , |3^ ], where FM states are found for doping regions with U/t>5. This 
discrepancy may be explained by the appearance of antiferromagnetic SRO which counteracts 
ferromagnetism (as in the SRO theories of finite-temperature magnetism, cf. Sec. 1) and which 
we have found to be stabilized for U/t > 5 — 6 (for details see below). The suppression of FM 
Hartree-Fock ground states by antiferromagnetic SRO is in accord with Quantum Monte Carlo 
(QMC) results ||. 

To discuss the competition between antiferromagnetic LRO and antiferromagnetic SRO 
and to give more insight into our SRO concept, already described qualitatively in Sec. 1, 
let us first quantify the behaviour of various magnetic parameters as functions of hole and 
electron doping. In Fig. 2 the sublattice magnetization itla in the AFM phase, the local- 
magnetization amplitude fh being finite only in the SRO-PM phase, where the SRO is found 
to be antiferromagnetic (Ji < 0), and the local magnetic moment m^oc = f ('^ ~ are 
depicted. In the case 6 < we have taken U/t = 9 (instead of U/t = 8 for 5 > 0), in order to 
get an SRO-PM phase at t' /t = —0.16. As can be seen, at large enough interaction strengths 
we obtain an AFM^SRO-PM phase transition of first order at a small critical doping 6c_^ 
(the thin dashed lines in Fig. 2 indicate in the region, where the AFM solution becomes 
metastable) and a SRO-PM^PM transition of second order at 6c^. 

The local magnetic moment measures the 'localization' of the electron spins and, contrary 
to m, is finite in the PM phase as well. With decreasing \6\ the correlation strength increases 
resulting in the increase of m/oc (as known from the theory of itinerant magnetism). If the 
'localization' prevails over the itinerancy, the formation of antiferromagnetic SRO is energeti- 
cally favoured {\5\ < \5c^\). Once the SRO has been established, mioc is slightly enhanced with 
respect to its value in the PM state (cf. Fig. 2). 

In the ground-state phase diagram of the 2D t-t'-Hubbard model shown in Fig. 3 one 
observes the suppression of LRO by SRO (at large enough U /t values) in a wide doping region. 
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Fig. 2: Local-magnetization amplitude rh (solid curve) compared with the sublattice magneti- 
zation rriA (dashed line) and the local magnetic moment rriioc (chain-dashed line) as functions 
of electron doping (U/t = 9; a,c) and hole doping (U/t = 8; b, d). Note the electron-hole 
symmetry for t' = and the asymmetry for t' ^ 0. The thin lines indicate the metastable 
solutions. 

The energy comparison between the SRO-PM phase and various incommensurate spiral LRO 
phases (obtained by the spin-rotation-invariant SB saddle-point solution [B3]) was performed 



in our previous work [12| for the particular value U/t = 8, where also the problem of phase 
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5 

Fig. 3: Ground-state phase-diagram of the t-t'-Hubbard model, where AFM, PM, and SRO- 
PM denote the antiferromagnetic LRO phase, the paraphase without SRO, and the paraphase 
with SRO, respectively. The solid (dashed) lines mark the phase boundaries at t' /t = —0.16 
{t'/t = 0). 

separation was discussed, and yields qualitatively the same SRO effects. For U/t = 8 (being 
realistic for LSCO) and S > we get 5c^ ~ 0.04 which agrees with the observed critical hole 
doping for the destruction of antiferromagnetism in LSCO. From this result we conclude that 
the rapid disappearance of antiferromagnetism in the quasi-2D hole-doped cuprates is related 
to the persistence of a strong 2D antiferromagnetic SRO in the paraphase. 

For t'/t = —0.16 and Q<U/t<8 the pronounced electron-hole asymmetry manifests itself 
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in the stability of the AFM phase in a much wider electron doping region as compared with 
the hole-doped case. The t'-induced shrinking of the stability region of the SRO-PM phase in 
the (5 < case originates from band-structure effects and, in particular, from the shift of the 
logarithmic Van-Hove singularity in the noninteracting density of states. Our findings are in 
accord with the asymmetry of the experimental phase diagram showing the disappearance of 
antiferromagnetism in Nd2-5Ce5Cu04 at 5c^ ~ —0.15. 

In the hole-doped systems, the hopping along the square lattice diagonals favours anti- 
ferromagnetic correlations (Fig. 3) so that the SRO persists up to 5c^ ~ 30% {U/t ~ 8). Since 
just at this doping the superconductivity in LSCO disappears, we suggest that the antiferro- 
magnetic SRO may be closely connected with the superconducting pairing mechanism in the 
hole-doped cuprates. 

Let us now discuss the weak-interaction part of the phase diagram of the t-f'-Hubbard 
model, where only the PM^AFM phase transition may take place. At half-filling and t' ^ 0, 
this transition occurs at a non-zero critical value Uc rather than at C/ = as for t' = 0; 
at t' It = -0.16 (-0.2) we get Ucji = 2.3 (2.5). This result agrees with the QMC value 
{Uc/t = 2.5 ± 0.25) obtained at t'/t = -0.2 [||, H and may be explained by the shift of the 
position of the Van-Hove singularity for t' ^ (for comparison, the Hartree-Fock calculation 
yields Uc/t = 2.1 [^]). A more detailed investigation of the half-filled i-t'-Hubbard model 
will be published elsewhere |Q. 

At the hole doping 6 = 0.13 and t'/t = —0.16 (see the phase diagram. Fig. 3), the AFM 
state is stable against the FM state down to C/ = which agrees with the Hartree-Fock result 
of Ref. {6 = 0.12 at t'/t = —0.2) but contradicts the result obtained by Brenig [^]. Thus, 
the correct weak-interaction limit of our saddle-point solution is confirmed. 

Let us stress again that at large interaction strengths (U/t > 5) the appearance of SRO, 
in addition to the inclusion of an appreciable part of the correlations at the SB-PM saddle 
point, changes the phase diagram qualitatively as compared with any Hartree-Fock solution, 
in particular, the FM and AFM LRO states are suppressed in a wide doping region. 
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5.2 Spin susceptibility and short— range order 

Based on the solution of the self-consistency equations discussed above, in the PM and SRO- 
PM phases we have evaluated the doping dependence of the zero-temperature susceptibility 
x(0, 5) according to the theory presented in Sec. 4. 

Figure 4 shows our results for the electron- and hole-doped cases at large interaction 
strengths, where the SRO-PM phase is stabilized (cf. Fig. 3). In the PM phase, the Pauli 
susceptibility of the SB-renormalized tight -binding band reveals a pronounced increase with 
decreasing (electron or hole) doping which notably depends on t' /t. This behaviour may be 
ascribed to Fermi-surface topology effects. Beyond the PM^SRO-PM transition at 6c^ (cf. 
Fig. 3), the Pauli susceptibility is suppressed due to the antiferromagnetic-SRO-induced spin 
stiffness against the orientation of the local magnetizations along the direction of the homoge- 
neous external magnetic field. Correspondingly, at 5c^ a cusp in x(0, 6) appears. According to 
the phase diagram, with increasing U /t the PM^SRO-PM transition shifts to higher doping 
levels so that the peak position in x(Oi'^) reveals the same U/t dependence. Analogously, from 
Fig. 3 the t' dependence of the peak position emerges. With decreasing doping the degree of 
SRO or, equivalently, the magnitude of the antiferromagnetic exchange energy, \Ji\, increases 



1 12 1 so that the susceptibility decreases. 

In the low-doping regime {\5\ ^ |(5c^|), where the SRO is most pronounced and the 'local' 
contribution to the susceptibility ( [4.1| ) predominates the 'itinerant' contribution, the absolute 
value of x(0, (5) may appreciably depend on the chosen approach in describing the local prop- 
erties for the following reason. In our scalar SB approach the local properties are treated on 
the basis of an effective Ising functional, describing longitudinal magnetization fluctuations 
only, which yields an Ising-type 'local' contribution to x- However, an improved SRO theory 
based on the spin-rotation-invariant SB scheme |2^ results in an effective Heisenberg func- 
tional including the transverse fluctuations too p7| ]. As known from the pure Heisenberg model 
p8| , the finite value of the zero-temperature susceptibility (being proportional to the inverse 
exchange coupling) is due to the transverse spin fluctuations. Therefore, we suggest that a 
spin-rotation-invariant theory of SRO describing the interrelation between a Heisenb erg-type 
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Fig. 4: Uniform static spin susceptibility at T = as a function of electron doping {U /t = 9; 
a) and hole doping {U/t = 8; b). 

'local' and 'itinerant' contribution to the susceptibility may enhance the magnitude of x in the 
low-doping region as compared with our results shown in Fig. 4. 

In Fig. 5 we have depicted the zero-temperature susceptibility calculated for the t-t'— 
Hubbard model taking realistic parameter sets for LSCO |17] without using any fit procedure 



in comparison with experiments on the spin contribution to the magnetic susceptibility of 
LSCO at 50 K |l^. This contribution is obtained from the experimental data on the 
total susceptibility ||l^ by subtracting the diamagnetic core (—9.9 x 10~^ emu/mol) and Van 
Vleck (2.4 X 10~^ emu/mol) contributions which, according to Ref. |16], can be taken as 
independent of doping and temperature over the limited parameter region studied here. As 
Fig. 5 shows, our SRO theory reproduces very well the position and, for suitable realistic 
parameters, even the height of the observed pronounced susceptibility maximum at a hole 
doping of about 25%. Moreover, the qualitative doping dependence of x is in reasonable 
agreement with experiments. However, in the low-doping region, the theoretical susceptibility 
is too low as compared with experiments. This deficiency may be ascribed to the neglect of 
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Fig. 5: Uniform static spin susceptibility at T = as a function of hole doping. The theoretical 
results obtained for the 2D Hubbard model at U/t = 6, t'/t = —0.16 (solid) and U/t = 8, 
t' /t = (chain-dashed), where the realistic value t = 0.3 eV is taken, are compared with 
the spin contribution (x) to the (corrected) experimental susceptibility of La2-5Sr5Cu04 at 
T = 50 K (Refs. [[| and Q). 

transverse magnetization fluctuations, as discussed above. 

Finally, we notice that the increase of the spin susceptibility upon doping, explained by 
our theory as an SRO effect in the one-band t-i'-Hubbard model at U/t > 5-6, is in qualitative 
accord with recent QMC data ||3^ and with results on the one-band Hubbard model {t' = 0) 



based on a semi-phenomenological approach |11] and on the composite-operator method 
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In those works a maximum in the spin susceptibihty was found even at a smaher coupUng 
{JJit = 4). 



6 Summary 

In this paper we have developed a theory of magnetic short-range order (SRO) in itiner- 
ant electron systems on the basis of the one-band t-t'-Hubbard model. The slave-boson 
functional-integral method was improved, as compared with previous approaches [1, 2, 4~6, 8, 
9] , by taking into consideration the self-consistent renormalization of the paramagnetic saddle 
point by SRO. The basic ingredients of our (non-CPA) theory are the following. 

(i) The four field slave-boson functional-integral representation is used to include an appre- 
ciable part of the correlations already at the uniform paramagnetic (PM) saddle point. 

(ii) The SRO of local magnetizations is incorporated beyond the PM saddle point within 
an effective generalized Ising model treated by means of the cluster variational method 
(CVM). The CVM is illustrated by cluster approximations with successively increasing 
sizes of the basic cluster. 

(iii) The SRO saddle point is determined from the free-energy functional by minimization 
with respect to all Bose fields and CVM variables. 

(iv) The inclusion of a homogeneous external magnetic field at the SRO saddle point allows 
the self-consistent calculation of SRO effects on the uniform static spin susceptibility. 

The theory is numerically evaluated at T = 0, where the SRO is taken into account at 
the Bethe-Peierls level of approximation. The dependence of SRO effects on the interac- 
tion strength, the hole and electron doping and on the band structure (modeled by t' /t) is 
investigated. The main results are summarized as follows. 

(i) At 5 < C//t < 10 {t' /t = —0.16) and for any doping level, the ferromagnetic long-range 



order (LRO) obtained in previous Hartree-Fock calculations [32] is suppressed due to the 
appearance of antiferromagnetic SRO. 
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(ii) The antiferromagnetic and incommensurate spiral LRO phases obtained by mean-field- 
type approaches make way to a paraphase with antiferromagnetic SRO in a wide doping 
region. In the ground-state phase diagram, a pronounced t'-induced electron-hole asym- 
metry is observed. 

(iii) The uniform static spin susceptibility increases upon doping in the SRO-PM phase and 
shows a cusp at the transition to the Pauli paraphase. 

From the good quantitative agreement of our theory with experiments on cuprates we conclude 
that the concept of a strong SRO in electron correlation models for high-Tc compounds may 
play an important role in the explanation of many unconventional magnetic properties and of 
the superconducting pairing mechanism in high T^s. 
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Appendix A: Effective Ising model 



According to the reasonings in Sec. 2.2, by the use of ( 2.1^ ) to (2.20) the functional ( p. 13 ) is 



expanded as 



(A.l) 



with 




(A.2) 



and the single-site and A'^-site fluctuation contributions 



'i^^'\{bs^,s,}) = j2{ud:d,^ 



(A.3) 



<J 



and 



'f'^^H{^,sJ) = 1 d^/(a;-/x)ilmTV In - Gl.^J^,%^,s^ 



(A.4) 
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respectively, where 



Using the identity 



we get 



i i'ij) (ijk) 



with 



hiiba,) = 



{"J 



Jij ipa - 5 ba ) 



Treating (A.4) in the pair approximation we have 



Tr In 



^ii-G'(^j)aTjAba^,aj 



^ In [l - G1ij^^Tj^{ba^,aj) GljiyT^„{ha^,ai 



(A.5) 



(A.6) 



(A.7) 



(A.8) 
(A.9) 
(A.IO) 



(A.ll) 



and the higher order terms in ( [A.7| ) vanish. Then, the functional ( |A.7] ) takes the form of the 
generalized Ising model ( 2.23D , where the effective parameters are given by 
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a 

)]}, (A.12) 

a n 

h{ba) = -^^a^Ud^d^ - n^iy^ + rh^^^ + ^(^0aa + l^Zn'Pn,aaa-) , (A.13) 
a an 

JniK) = -\'^{'^n,aaa -^n~aaa) ■ (A.14) 

Here, <P^„ = ^i^(aj|^^^^, and 'Pn,aa'a = <P(ij)^^{ai, aj) ,^=<, with 

a , —a' 

^fe)„- = / rf'^/(^-/^)^Im ln[l-q,^.)^,r,.,(6„^.,a,0%)„.r,J (A.15) 
is the two-site fluctuation contribution which takes into account SRO effects. 

Appendix B: Cluster variational method and examples 

To ihustrate the algorithm of the CVM sketched in Sec. 2.3, we consider the simplest cluster 

approximations with the basic clusters (a), (b) and (c) of Fig. 1 corresponding to the Bragg- 
Williams |41|, Bethe-Peierls [20| and Kramers-Wannier approximations, respectively. The 



(v) (v) 

energies Ek and their multiplicities Ak for the topologically inequivalent spin configurations 



of the isolated clusters are given in Table 1. The consistency relations ( 2.27] ) between the distri- 



bution numbers pl^^ of the subclusters and the basic clusters p^T^ as well as the dependences of 
pli'^ on the mutually independent CVM variables for the (t(;)-clusters are listed in Table 2. 
Concerning the choice of the variables in proceeding from the (w)-cluster to the {w + 1)- 
cluster approximation (e.g. (a) (6)), it is convenient to preserve the definitions of the 
or, equivalently, the x-dependences of plT'* with (v) C [w). Then, only the x~dependences of 
P^k'^^^ have to be determined from the consistency relations. Accordingly, we take the choice 

(a) (a) (b) (c) (c) , (c) (c) 

As can be seen from Table 2, the expectation value Kj^ = (sj) is given by Kj^ = Xi (being 
the order parameter in the ferromagnetic Ising model). 

The geometrical factors yio ^appearing in eq. ( |2.24| ) are defined by the set of equations |21] 



^a(")('')yl'^)=x("); {u)^{w), (B.l) 
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Tab. 1: The topologically different configurations k (• =|, o =|), multiplicities aI^^ and 
configuration energies el^^ for the clusters (a), (6) and (c). 

where the matrix element a^'^^^'"^ indicates the number of clusters (n) contained in the cluster 
(v) and x^"^^ gives the total number of (u)-clusters in the lattice divided by A^. 

Let us now exemplify the CVM in the course of successively increasing sizes of the basic 
cluster (w). 

(i) Bragg-Williams approximation 

Taking a single site as basic 'cluster', i.e. {w) = (a), from Tables 1 and 2 we obtain the 
free energy ( 2.25| ) as 

2 



^^(xi) = -iV 



' K = l 



(B.2) 



Since there are no subclusters and the total number of (a)-clusters in the lattice coincides 
with N = 1), from ( |B.1| ) we get trivially yi'^^ = 1. Thus, the free energy ( |2.24 ) in 
the Bragg-Williams approximation is 



^„(xi)=^(")(xi). 



(B.3) 



The CVM equations (|2.28|) reduce to the equation jj^^ = 0, which, in the paraphase 
and the h = limit, has the solution = 0. 
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(c) 


p'q'^ = i(l-xi-4x2 + 2x3 + 2x5) 


(b) 




pS'') = i(l + Xi-2x2) 




(^') (c) , (c) , (c) , (c) 


P2 =^2 






Pf = W-^l- '^^2) 


(a) 


(a) (b) , (b) 
Pi =Pl +P2 


p('^) = l(l + xi) 




(a) (6) , (fe) 


P2^ = W-Xi) 



Tab. 2: Distribution numbers p^f^ of the configurations k (Table 1) as functions of the distri- 
bution numbers of the basic cluster [w) (consistency relations) and as functions of the CVM 
variables x^. Note the hierarchy in the choice of x^ (see text). 



X 



,h + {l-Ax2)J,-Wp^^^^^p''^ 



(B.4) 



(ii) Bethe-Peierls approximation 

Choosing the cluster (6) as basic cluster (Table 1), the number of relevant spin configu- 
rations and CVM variables (Table 2) increases by one. Then, the free energy (I^H) of 

an assembly of N independent Bethe clusters (6) is obtained as 

3 

Determining the geometrical factors fjw^ from (|B.1| ), we note that the cluster (a) is con- 
tained twice in (6), and the cluster (6) is contained in the lattice 2N times. Accordingly, 
the equations (lB.l| ) become 



(B.5) 



K = l 




28 



From the solution of ( |B.5D and by ( |2.24| ) we obtain the free energy in the Bethe-Peierls 
approximation, 

J6(xi,X2) = 2^W -S^^'^) , (B.6) 
with j^(^) and given by ( [B.4[) and (|B.2| ), respectively. 

By ( p. 31 ), the nearest-neighbour pair correlation function K j = {s^Sj)i is expressed as 



Kj = 1-4x2 . 



Concerning the CVM equations j^^^ = {l = 1,2), the equation for xi reads 

1 



(B.7) 



h + 



3 , 1 - xi , 1 + xi - 2x2 

— m h m 

2 1 + xi l-xi-2x2 



. 



(B.8) 



In the paraphase, eq. ( [B.8D has a non-trivial solution for h ^ only and therefore is 
needed in the calculation of the spin susceptibility (Sec. 4). In the h = limit (xi = 0) 
there remains the second CVM equation 



with 



The solution gives 



2x2/f = 1-2x2 



/„ := exp{pjn} 



X2 = i[l — tanh {pJi 



(B.9) 



(B.IO) 



(B.ll) 



Inserting (B.ll) into ( |B.7] ) one gets the same correlation function Kj as in the original 
Bethe cluster method (working with an effective Bethe field). 

(iii) Kramers- Wannier approximation 

Here, the cluster (c) is taken as a basic cluster, and the number of CVM variables is 
enlarged by three (Tables 1, 2). The free energy (p.25|) for (v) = (c) is given by 



J^('^)(xi,...,X5) = -7v{4xi/i + 4(l-4x2)Ji 

+2[1 - 4(X2 + X3 - X4)] J2 - ^ E A(^)p(^) lnp(-)} . (B.12) 



K = l 



29 



From ( [B.l| ) we obtain the geometrical factors as solutions of the equations 



1 2 4 
14 
1 



yc 



yc 

V^^ 7 



1 

2 



(B.13) 



which yield the free energy in the Kramers-Wannier approximation, 



. . . , xs) = T^"^ - 2^W + ^('^) . 



(B.14) 



Besides the nearest-neighbour correlation function ( |B.7] ), from ( p.31| ), ( |B.12 ) and (B.14) 
the next-nearest-neighbour correlation function 



Kj^ = 1 - 4(X2 +X3- X4) 



(B.15) 



can be calculated. In the paraphase and in the h = limit, the CVM equations yield 
xi = X5 = 0, and the remaining equations become 

P^^P?l!=P^2''p? , P^s^ll=pi'\ P^MVf=pi'^' (B.16) 



with the distribution numbers taken from Table 2 and the definitions ( B.lOj ). Instead of 
giving the explicit solutions of ( [B.16| ), we sketch the resulting zero-temperature behaviour 
in Fig. 6. There, the degenerate ground-state spin configurations for h = are shown 
in the J1-J2 plane. For J2 = 0, we have antiferromagnetic (vr, vr) SRO if Ji < and 
ferromagnetic (0, 0) SRO if Ji > 0. The inclusion of the next-nearest-neighbour Ising 
exchange energy J2 gives rise to frustration effects and permits to describe SRO of a 
stripe structure denoted by (vr, 0)-SRO. 
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Fig. 6: Paramagnetic solutions of the CVM equations at T = within the Kramers-Wannier 
approximation. 
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